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ABSTRACT 

We investigate faint radio emission from low- to high-luminosity Active Galactic Nuclei (AGN) 
selected from the Sloan Digital Sky Survey (SDSS). Their radio properties are inferred by co- 
adding large ensembles of radio image cut-outs from the FIRST survey, as almost all of the 
sources are individually undetected. We correlate the median radio flux densities against a range 
of other sample properties, including median values for redshift, [O III] luminosity, emission 
line ratios, and the strength of the 4000A break. We detect a strong trend for sources that 
are actively undergoing star-formation to have excess radio emission beyond the ~ 10 28 ergs s _1 
Hz -1 level found for sources without any discernible star- formation. Furthermore, this additional 
radio emission correlates well with the strength of the 4000A break in the optical spectrum, and 
may be used to assess the age of the star-forming component. We examine two subsamples, one 
containing the systems with emission line ratios most like star-forming systems, and one with 
the sources that have characteristic AGN ratios. This division also separates the mechanism 
responsible for the radio emission (star- formation vs. AGN). For both cases we find a strong, 
almost identical, correlation between [O III] and radio luminosity, with the AGN sample extending 
toward lower, and the star-formation sample toward higher luminosities. A clearer separation 
between the two subsamples is seen as function of the central velocity dispersion a of the host 
galaxy. For systems at similar redshifts and values of <r, the star-formation subsample is brighter 
than the AGN in the radio by an order of magnitude. This underlines the notion that the radio 
emission in star-forming systems can dominate the emission associated with the AGN. 

Subject headings: galaxies: active — galaxies: starburst — radio continuum: galaxies 



1. Introduction Ikauffmann fc Haehneltl200oH iHeckman et al . 2004) 



There are several observational lines of evidence 
The formation of the bulge component of galax- for this theoretical claim. For instance, the cen- 

ies and their central black holes is now understood tral ve i ocity dispersion of bulges is found not just 

to be tightly linked (e.g., | Richstone et al. | [l998|; to correlate with the mass of the bulge, but also 
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with the inferred mass of the c e ntral black hole 
(e.g., Ferrarese &: Merritt 200d ; Gebhardt et al. 
2000). Since the correlation between the mass of 



the central black hole and the bulge exists over 
a large range of bulge masses, it is thought that, 
above some lower mass limit, all bulges harbor 
massive central black holes. Locally, accurate 
black hole mass measurements have been made 

200.4 



Schodel et al 



for our own G alaxy (e.g. 

iGhez et"aL 20051). and t he Andromeda galaxy 
(e.g.. iBender et all 120051 ). Our own Galactic 
black hole is rather dormant ba sed on its rel- 
ative l y low X-ray em ission (e.g., iBaganoff et al 



200.4 ku et aj]l200ri ). but in general these mas- 



sive black holes arc considered a key component 
of Active Galactic Nuclei (AGN), in which accre- 
tion onto the black hole itself can p roduce large 
amounts of energy efficiently (e.g., Lvnden-Belll 
19691 : IPringlelll98lh . 



It is of interest, then, to study AGN properties 
over a large range of luminosities, ranging from 
the AGN-emission-dominated quasars to galaxies 
for which the presence of a weak AGN can be in- 
ferred through unusual optical emission line ratios. 
Kauffmann et al. (2003a, b) studied the proper- 
ties of the host galaxies of AGN, using large, well- 
defined samples from the Sloan Digital Sky Sur- 
vey ( SPSS, e.g. JYork etHlbOOOHStoughton et al 



2002). The sheer size of these samples (a few tens 



of thousands of galaxies) allows for very accurate 
measurements of sample averages and their trends 
as function of optical continuum and emission-line 
(mainly [O I II]] luminosities. One of the con- 
clusions from iKauffmann et al. ( 2003aL hereafter 
K03a) is that the main difference between low- and 
high-luminosity AGN is the significant presence of 
a young stellar population in the latter. 

In this paper, we will concentrate on the radio 
emission properties of a large sample of AGN, sim- 
ila rly selected from the SPSS survey (Pata Release 



4. lAdelman-McCarthv et al.l l20odb . We use the 



latest AGN compil ation of Kauffmann et al[j] who 
used the BPT (jBaldwin. Phillips, fc Terlevicb 



19811 ) diagram of emission line ratios to select 



for AGN activity. The current version of this cat- 
alog contains 80,156 candidate objects which are 
also covered by the FIRST radio survey 



(jBecker. White, fc Helfandlll995l ). In addition, we 
selected 14,165 quasars from the PR4 release for 
which we have FIRST images and [O III] line 
coverage (limiting the quasar redshifts to < 0.8). 
These quasars will form the high-luminosity por- 
tion of our sample. 

We focus our attention on the radio emission 
of the AGN for several reasons. First, almost all 
AGN are thought to produce radio flux at some 
level; even the Galactic black hole, currently in a 
"quiescent" phase, emits in the radio. Some AGN 
are extremely luminous in the radio (compared to 
their optical output) - the so-called radio-loud ob- 
jects, although we are not considering these specif- 
ically, as our average AGN is radio-quiet by a wide 
margin. Secondly, radio emission is unaffected by 
intrinsic absorption within each source. It there- 
fore provides a more or less accurate measure of 
AGN radio output, with the sole caveat being rela- 
tivistic beaming effects. Since none of our sources 
are classified as Blazars for which beaming effects 
can be large, we assume this to play a minor role. 
Finally, the ready availability of a large-area, high- 
resolution radio survey in FIRST provides radio 
information on the large numbers of AGN in its 
9033 deg 2 coverage. 

There is, however, one apparent problem: al- 
most none of the sources in the SPSS sample are 
detected by the FIRST survey, with its flux den- 
sity threshold at 1.4GHz of 1.0 mJy (5a). Out 
of 80 154 AGN, only 5 875 (7.3%) are detected di- 
rectly by the FIRST survey. The relevant numbers 
for the quasars are only slightly better: 1 493 out 
of 14165 (10.5%). Clearly by limiting the sam- 
ple to the brightest sources, one is ignoring the 
bulk (~ 90%) of the population. Fortunately, we 
can stack radio images, and detect the mean and 
median peak flux densities of ensembles of unde- 
te cted sources. This m ethod is discussed in depth 
bv I White et all (|2006l . hereafter Paper I); we sum- 
marize the method in § [21 Given our large sample 
size, we are able to detect not just the (median) en- 
semble flux density values, but also discern small 
trends in those values. For our typical stack sam- 
ple size of 5000 sources, we attain an rms noise of 
2.6/xJy, and are therefore very sensitive to small 



1 http: / / www.mpa-garching.mpg.de/SDSS /DR4 

2 The Very Large Array is an instrument of the National 



Radio Astronomy Observatory, a facility of the National 
Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 
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changes in the sample median values as we vary 
other parameters such as the redshift distribution 
and the [0 III] luminosity distribution. We only 
have a single sample, so the way we investigate the 
varying dependencies is by re-grouping the stacks 
according to the quantity of interest. For instance, 
if we order the sample by redshift and generate 15 
adjacent bins (5343 sources per bin), the result- 
ing radio stacks will be different from the 15 bins 
sorted by [O III] luminosity; each plot in this pa- 
per contains the exact same data, but arranged 
differently. 

The typical median radio flux densities for these 
sources is found to be on the order of 50 to 
100 /uJy, well below the FIRST threshold, but 
clearly within the capabilities of the VLA. Indeed, 
various very deep, small-scale radio surveys ex- 
ist (e.g., Hopkins et al.lll998t iRichards et al.|[l999l 



de Vries et al.l 120021) . At these low flux densities, 



star-formation could account for a large fraction of 
the radio emission in our objects, since flux density 
levels of 50 /iJy are well within range of starform- 
ing galaxi es (without AGN) at s imilar redshifts 



see, e.g., iWindhorst et al.l 11999c iFomalont et al 



2002). We cannot, therefore, assign all of the ra- 
dio emission to the AGN, as some fraction of the 
emission may be due to star-formation. 

The paper is organized as follows. In § [3 we 
summarize the stacking method as presented in 
Paper I. The next section describes the proper- 
ties of the complete sample of 80,156 sources. We 
then use emission line ratios as well as the strength 
of the 4000A break to isolate two subsamples of 
26,715 sources each (§ [3]). The hrst subsamplc 
has properties most consistent with the presence 
of an AGN (the "pure" AGN sample); the second 
has star-formation characteristics while still meet- 
ing the AGN selection criterion of K03a. Section 
2] details the differences and similarities of these 
subsamples, which are discussed in further in § O 

2. FIRST image Stacking Technique 

Paper I demonstrates that it is possible to mea- 
sure the mean and median radio flux density val- 
ues of distributions of sources, even though indi- 
vidual sources fall far below the detection thresh- 
old of the FIRST survey. Provided that the sam- 
ple is large enough, one can attain rms noise values 
in the radio sky well below the canonical FIRST 



value of 0.15 mJy. Actual snap-shot stacking ex- 
periments of blank pieces of skj{f| conform to the 
expected l/VN behavior in the pixel statistics, 
down to better than 1/iJy. Based on these num- 
bers, it is clear that one can detect stacked point- 
source mean flux densities of a few tens of /iJy 
with high fidelity. 

As outlined in Paper I, we prefer to use the me- 
dian value of the distributions over the mean. The 
latter quantity is rather easily affected by outliers 
with large flux densities (the ones which are actu- 
ally above the FIRST detection threshold). Sim- 
ply removing the sources above the threshold from 
the sample to arrive at the mean of the undetected 
sources is not robust, as a small change in the cut- 
off flux density (e.g., from 1.0 to 0.9 mJy) results 
in a significant change in the mean flux value. 

The second concern addressed in Paper I is the 
calibration of the stacking procedure for effects in- 
troduced by the radio data analysis in general and 
the 'CLEAN' algorithm in particular. Analysis 
of both stacked artificial sources, and actual sub- 
threshold sources detected in full-synthesis, deeper 
radio imaging shows that one does not recover 
fully the flux that went into the stack. These re- 
sults are illustrated in Fig. 2 of Paper I; the cor- 
rection for this "snapshot-bias" for sub-threshold 
sources is given by S P}Corr — 1.4QS p , where S p is 
the median peak flux density value. All of our 
stacked peak flux density measurements have been 
corrected by this factor. 

2.1. Stacking luminosity images 

Since most of our subsequent discussion deals 
with radio luminosities, and not radio flux den- 
sities, we have to consider whether there are any 
distribution peculiarities that might affect these 
quantities. We do not know the flux density dis- 
tribution of sub-threshold sources, nor the na- 
ture of the sources that make up this popula- 
tion. Below flux densities of a few mJy, the com- 
position of the radio source population changes 
from AGN -dominated to star-fo rmation domi- 
nated (e.g.. lWindhorst et al.lll999h . 



3 No pieces of the sky are truly blank, but will consist of faint, 
unresolved background objects. Our stacking technique is 
not sensitive enough to pick up this signal: we measured 
a background signal consistent with 0/xJy to within the 
0.9^tJy rms noise (for 80 154 empty patches). 
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There are two ways of calculating the median 
radio luminosity. One is simply to stack the cut- 
outs (with pixels in units of Jy), and use the me- 
dian peak flux density and median redshift (of the 
sources that were used in the stack) to calculate 
the median radio luminosity. The tacit assump- 
tion here is that those two median values actually 
describe the same "median object" . The alterna- 
tive approach is to convert each snapshot into a 
luminosity image using the redshift of the individ- 
ual source, and then stack these images instead. 
Each pixel in a snapshot image is converted from 
flux density to luminosity using the following cos- 
mological parameters: H Q = 70 km s _1 Mpc -1 , 
= 0.7, and Qm = 0.3. 

It turns out that the results are not exactly the 
same. Table 2 (columns 4 and 5) lists the radio 
luminosities derived using both methods for the 
sample sorted in D n (4000) strength, showing typ- 
ical differences of less than 25%. For the sake of 
consistency, we will use the luminosity-stacked val- 
ues instead of the flux density-stacked ones for the 
remainder of the paper. 

3. Star-formation vs. AGN activity 

Our data sample was initially selected based on 
source locations in the BPT plot (Kauffmann et 
al. 2003a, see also Kewley et al. 2006). The ra- 
tios of the narrow-lines of [N II] and [O III] over 
the permitted Balmer lines Ha and H(3, respec- 
tively, serve as an excellent proxy for the amount 
of st ar-formation in an object (e.g., lOsterbrock 
1989). Our galaxies with AGN cover a large range 
in the BPT plot, and are thus affected by star- 
formation to varying degrees; the objects found 
close to the sample cut-off arc, toward the left of 
the distribution (see Fig. [1]) , have optical spectra 
dominated by emission line ratios commonly asso- 
ciated with ongoing star-formation. It should also 
be noted that Kauffmann et al. used a more liberal 
cut for their AGN c lassification than for instance, 
Kewlev et al.l (|200ll ). resulting in a fraction of the 
sources being more H II region-like than AGN-like. 
These sources are classifie d as "c omposite" in the 

( 20041 ) paper. 



related Brinchmann et al 



We use the 4000A break strength inde x D n (4000) 
as defined by I Kauffmann et al. ( 2003b ). This in- 
dex can be parameterized as the age since the 
last (instantaneous) episode of star-formation 




-0.5 o 
log([NII] / Ha) 



Fig. 1. — Emission line diagnostic plot (BPT, 
Baldwin. Phillips, fc Terlevichl 



19811 ) which dis- 



tinguishes between star-formation-dominated and 
AGN-dominated systems. The sharp cutoff re- 



flects t he selection boundary of Kauffmann et al 



(|2003al ). Points to the left of this cutoff are con- 
sidered pure star-formation without an AGN con- 
tribution. Pure AGN ratios are only expected in 
the extreme upper right hand corner, so most ob- 
jects are mixes. The grey scale squares illustrate 
the progression of low-strength D n (4000) systems 
(i.e., most actively star forming) to large break 
systems (light to dark grey respectively). The cir- 
cled groups, I and II, are defined in Fig. O 



(see Fig. 2 of Kauffmann et al. 2003b), in the 
sense that smaller break strengths correspond 
to more recent star-formatiorjj]. The ranges of 
the D n (4000) under consideration for our sources 
roughly correspond to ages since star-formation 



between 10 s 



6-10.0 



years. The short end of this 



range is comp arable to typical galaxy merger 
times cales Ce.g.. lSpringel. Di Matteo. fc Hernquist 
2005), whereas ~ 10 10 years at z ~ 0.2 is close to 



the age of the universe and therefore represents the 
initial burst of star-formation. Since we are not go- 



4 It should be noted that Kauffmann et al. (2003a) found 
that the AGN contribution to the blue continuum, which 
would give rise to an overestimate of the star-formation 
rate using this method, is small enough to be ignored for 
these objects. 
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Fig. 2. — Radio Luminosity as function of 
D n (4000) strength. The bins, each containing 
about 5000 sources, are sorted in order of increas- 
ing D n (4000) strength (from left to right), and are 
coded for [O III] luminosity. The D n index is not 
very dependent on redshift; the median redshifts 
range from 0.12 down to 0.07 (from left to right 
in the plot) . The age indicated on the top of the 
plot is the inferred age since the last burst of star- 
formation (in years), based on a linear conversion 
between D n (4000) and age (see text). The sub- 
samples marked I and II are described in the text. 



Fig. 3. — Correlation between fiber aperture- 
corrected star-formation rat e, based on the 
stren gth of the D n (4000) break IjBrinchmann et al 



20041 ) , and the radio luminosity. The circles repre- 
sent our sample, and the small triangles are star- 
forming galaxies without an AGN (taken from 
Brinchmann et al.). The least-squares fit to the 
non-AGN data points (black line) is given by: 
log(L£, on - AGN ) = (1.37±0.02)log(SFR) + (27.67± 
0.01). The red line represents the sum of the 
non-AGN component £p° n ~ AGN arLC [ a constant 
4.5 x 10 27 ergs s _1 Hz^ 1 due to the AGN. 



ing to distinguish between starbursts of different 
metallicities, we are approximating the relation 
between age and D n (4000) by a straight line over 
the D n (4000) range 1.3 to 2.0 (the relevant range 
for our sample): log (age) = 6.21 + 1.87D n (4000), 
with age given in years. 

The results for median radio luminosity as func- 
tion of the strength of the 4000A break are plot- 
ted in Fig. [2] It is clear that the brightest ra- 
dio emission is associated with the smallest breaks 
(i.e., the most recent episodes of star-formation). 
The approximate time since the last epoch of 
star- formation is given across the top of the plot. 
Above a D n (4000) index of about 1.8, there is 
no change in radio luminosity, which levels off at 
about 10 28 ergs s _1 Hz -1 . This might reflect the 
absence of any contribution to the radio emission 



from star-formation. The color coding indicates 
the [O III] luminosity, which is also seen to decline 
steeply as D n (4000) increases. 

3.1. Line emission diagnostics 

As a double check on the level of star-formation, 
we can see where the D n (4000)-sorted data points 
from Fig. [2] land on the BPT plot. In other words, 
do the sources with the low D n (4000) values have 
emission line ratios that reflect a higher incidence 
of star- formation? The results are plotted in Fig.Q] 
as a grayscale sequence of squares, running from 
light-gray to black for the small-to-large valued 
D n (4000) systems. It is clear that these two quan- 
tities do correlate rather well for our objects, and 
trace out a sequence of decreasing star-formation 
versus increasing 4000A break strength. The gray 
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squares in the top right corner have emission line 
ratios completely consistent with pure AGN emis- 
sion, whereas the squares to the left are signifi- 
cantly affected by star-formation. 

The combination of Figs. [1] and [2] leads us to 
suggest that for the "pure" AGN cases as indi- 
cated by large 4000 A breaks, the radio luminos- 
ity is strictly due to the AGN with average lumi- 
nosities near fO 28 ergs s _1 Hz -1 . The additional 
radio luminosity we see for other sources can be 
attributed completely to ongoing star-formation. 
This star-forming component can be an order of 
magnitude more luminous than the radio compo- 
nent due to the AGN. 

iBrinchmann et aT (2004) use the strength of 
the 4000A break to calculate the inferred star- 
formation rate (their sample not only encompasses 
our AGN sample, but also non-AGN which land 
to the left of the distribution in Fig. QJ. They cor- 
rected for the limited angular extent of the SDSS 
fiber which typically does not contain all of the 
light of the galaxy. We use these corrected val- 
ues to sort our sample in ascending order of star- 
formation rate. The median stacked radio lumi- 
nosities for 15 bins of ~ 5000 sources each are cal- 
culated and plotted against median star-formation 
rate in Fig. [3l The strong correlation again under- 
lines the importance of the star-formation com- 
ponent to the overall radio luminosity, even for 
systems containing AGN. The 10 28 ergs s _1 Hz -1 
"ground-state" radio luminosity from Fig. [5] cor- 
responds, based on this plot, to a star-formation 
rate of ~ 1 M yr _1 , or less, whereas the largest 
radio luminosities from Fig.[2]have rates exceeding 
10 Mq yr" 1 . 

The non-AGN portion of the Brinchmann et al. 
sample, once sorted by star-formation rate, exhibit 
a very tight correlation with median radio lumi- 
nosity (Fig. [3|). The least squares fit is given by 
log(L£° n - AGN ) = (1.37±0.02)log(SFR) + (27.67± 
0. 01). This sl ope of 1.37 is almost identical to that 
of lBelil (|2003l ) who find a slope of 1.30 for radio lu- 
minosities less than 6.4 x 10 28 ergs s _1 Hz -1 . The 
Bell results, however, show a different normaliza- 
tion, suggesting they underestimate the amount 
of star-formation for a given radio luminosity by 
a factor of 2. 

Even more remarkable than the tight correla- 
tion between estimated star-formation rate and 
mean radio luminosity over two orders of mag- 




Redshift 

Fig. 4. — Distribution densities for the 5 lowest 
(red) and 5 highest (blue) radio luminosity bins of 
Fig. [2] While the redshift coverage is comparable, 
there is a clear dependency on the strength of the 
D n (4000) index in the sense that it is easier for a 
given burst of star-formation to "fill in" the 4000A 
break in a smaller, low a system. The density 
contours are given by 2", with n = 2, 3, 4, • • •. 



nitude is that the AGN contribution to the ra- 
dio luminosity can be represented by a constant 
(4.5 x 10 27 ergs s^ 1 Hz -1 ) which, once added to 
the star-formation estimate (black line), gives us 
the red line. This single constant implies that the 
AGN output (in the radio) is more or less indepen- 
dent of the star-formation level of the host galaxy. 
It is also a better measure of this AGN "ground- 
state" than the 10 28 value derived from Fig. [5] 

In summary, since the star-formation rates are 
based on the D n (4000) values, we will use the lat- 
ter to isolate two subsamples; one that presumably 
has no, or very little, ongoing star-formation (la- 
beled I in Figs. [T]and [2]), and one that contains 
the systems with the most ongoing star-formation 
(labeled II) . Each one of these subsamples contain 
26,715 sources. 
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4. A closer look at the "pure" AGN sub- 
sample 

We designate the five rightmost data points in 
Figs Q] and [2] as representing the most AGN-likc 
sources. Their emission line ratios are the least 
star-formation-likc, and their median [O III] and 
[O II] luminosities are among the lowest in our 
sample (see Table 2). 

Figure [4] shows the distribution densities in red- 
shift versus velocity dispersion space for the large 
4000A break, low star-formation sources (red con- 
tours), compared to the sources with the small- 
est 4000A break and the highest levels of star- 
formation (blue contours). The redshift distribu- 
tion for each subset is comparable, whereas there 
is a clear offset in velocity dispersion. This is most 
likely due to the way the subsets have been se- 
lected. For a given burst of star-formation with 
some upper limit to its optical luminosity, it will 
appear more prominent (with a consequently lower 
value of the D n (4000) index) in smaller galax- 
ies than in large systems. Since the mass and 
size of the galaxy spheroi d scales with its cen- 
tral velocity dispersion (e. g. , iGebhardt et al 
Ferrarese fc Merritt 200(1) . we will be biased to- 
ward finding more star-formation among the lower 
a systems. 

4.1. Radio Luminosity and redshift depen- 
dence 

Luminosity, as in all flux-limited samples, is 
closely correlated with redshift, resulting in an ap- 
parent correlation that is dominated by the lack of 
bright sources locally, and faint sources at higher 
redshifts (due to the detection limit). This is il- 
lustrated in Fig. [5] Both subsamples I and II were 
sorted in redshift, and divided up into 13 bins of 
2056 sources each. The "pure" AGN-like bins are 
plotted as pentagons, and the sources most af- 
fected by star-formation as stars. In addition, we 
plotted 7 bins of 2025 quasars each, depicted as 
green squares. 

It is immediately clear that the median radio 
luminosities of the quasars are directly compara- 
ble to the narrow-line AGN values. There is no 
a-priori reason for the quasar radio luminosity me- 
dians not to be much brighter than the galaxies, 
but their values are comparable at the same red- 
shift. The quasar radio luminosities are slightly 



SF dominated 
£ AGN dominated 
■ Quasars 




Redshift 

Fig. 5. — Radio luminosity as function of red- 
shift. The narrow-line AGN sample is separated 
into AGN dominated (I, pentagons) and star- 
formation dominated (II, stars) parts. In addi- 
tion, 7 low-redshift quasar bins (squares) are over- 
plotted. The least squares fits are: log(Lri) = 
(1.98±0.07) log (z) + (30.75±0.07), and log(L R ) = 
(2.45 ±0.11) log \z) + (30.81 ±0.12) for subsamples 
I and II respectively. The red solid line indicates 
the radio luminosity of a constant 60 fi3y source. 

less than the star-formation values, and slightly 
more than the AGN subsample values (at the same 
redshift). On the other hand, quasars are con- 
sistently brighter in [O III] luminosity than the 
narrow-line AGN at similar redshifts (see Fig. O . 

The fact that the median radio luminosity is 
the same would suggest that we are observing an 
isotropic emission component, and not orienta- 
tion / obscuration-dependent quantities like emis- 
sion lines and the presence of broad permitted 
lines. The median radio luminosity therefore un- 
derlines the vali dity of the standard AGN unific a- 
tion model (e.g.. lBarthelll989t lAntonuccilll993h . 

Another result based on Fig. [5] is that the star- 
formation-dominated galaxies are brighter in the 
radio than the pure AGN galaxies, albeit that the 
difference becomes smaller at redshifts beyond 0.1 
(where both samples I and II merge into the low 
end of the quasar trend). The more or less con- 
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[0 III] Luminosity log(ergs/s) 
39 39.5 40 40.5 41 41.5 42 42.5 




[0 III] Luminosity log(L/L s ) 



Fig. 6. — Correlation of median values for radio 
luminosity and [O III] luminosity. The color in- 
dicates the median redshift for a particular bin. 
The 7 squares are broad-line AGN / Quasars. The 
units on the top and right hand side are converted 
from the measured values by the following con- 
stants: log (1.4 x 10 9 ) Hz, and log (3.826 x 10 33 ) 
ergs/s for the x-axis (observing frequency) and 
y-axis (bolometric luminosity Sun) respectively. 
Note that for all these points the radio luminosity 
is at least 3 orders of magnitude lower than typ- 
ical radio-loud objects. Symbols are the same as 
in Fig. 03 

stant offset of a factor ~ 3 (0.5 in log) between the 
two subsamples (below z ~ 0.1) suggests that the 
average radio luminosity due to the AGN accounts 
for at most 25% of the total radio emission during 
an ongoing phase of star-formation. 

4.2. [O III] luminosity as AGN activity in- 
dicator 

The [O III] emission-line luminosit y offers an- 
other w ay of gauging the AGN output. 
(<2003al) estimate that perhaps as little as 7% of 
the [O III] line emission can be attributed to star- 
formation in their AGN composite spectrum. This 
means that, even though both the radio and [O III] 
luminosities depend on redshift, their correlation 
should be relatively free of redshift bias (by drop- 



ping out of the ratio) . 

For this purpose, we first order the sample in 
[O III] luminosity, and then stack the cutouts to 
calculate the median radio luminosity The [O III] 
luminosities have been corrected for extinction us- 
ing the Ha / H/3 emission line ratio, under the 
assumption that the [O III] emission is co-spatial 
(see K03a). The [O III] sorting and stacking also 
helps in minimizing the redshift differences among 
the bins. Ideally, one would like to have similar 
source redshift distributions for each [O III] lumi- 
nosity bin, which would in effect isolate one from 
the other. 

The results are presented in Fig. [SJ We see a 
strong correlation between [O III] and radio lu- 
minosity. For comparison, we also included seven 
bins of low-redshift quasars. The spectral coverage 
of the SDSS is such that the [O III] emission line is 
included in the spectrum up to redshifts of ~ 0.8. 
This results in a sample of 14,165 DR4 quasars 
that have z < 0.8 and are covered by the FIRST 
survey. We split the sample into seven equal-sized 
bins, sorted in [O III] luminosity. Since we are 
looking at presumably unobscured quasars, we do 
not correct the [O III] luminosities for intrinsic ex- 
tinctior|f|. 

The two solid lines are least squares fits to the 
quasars and the star-formation subsample. The 
first ten data points of the AGN subsample (pen- 
tagons) are more or less consistent with the quasar 
relation, though we do note that there are excur- 
sions along the quasar fit. This may have to do 
with the extinction correction being less than ideal 
for the AGN subsample. The star-formation sub- 
sample does not exhibit similar behavior. The fits 
are: 

ir.AGN = (0.87 ± 0.02)L [O iii] + (2.5 ± 0.8) 
L r , SF = (0.70 ± 0.03)L [ora] + (9.1 ± 1.2) (1) 

with all the luminosities in units of log(ergs/s). 
It should be noted that the origin of the radio 

'We know that the main sample of galaxies has to have some 
extinction towards the AGN (they would otherwise appear 
as type 1 Seyfert objects), so correcting for the [O III] 
luminosity (mainly due to the AGN) makes sense. But 
since the quasars are type 1 objects, we will assume that the 
corrections for the quasars are considerably smaller than 
the typical factor ~ 10 corrections in the median [O III] 
luminosity for the galaxies. 
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emission is different for the two subsamples. For 
the AGN- dominated sample the radio emission is 
largely due to the AGN (though the star-formation 
component dominates at the highest luminosities, 
see Fig. [3]), whereas the AGN component is only a 
small contributor for the star-formation subsam- 
ple. As such, it is clear that the differences in the 
[O III] - radio luminosity relations cannot be ex- 
plained by a simple translation along either the 
[O III] or the radio luminosity axes. From Fig. [5] 
we already know that the starforming subsample 
is brighter in the radio than either the AGN or 
quasar samples (at the same rcdshifts). Therefore, 
the correct conclusions from Fig. [S] are: the star- 
forming subsample is both brighter in [O III] and 
in radio luminosity than the AGN subsample at 
the same redshift, and they have a different slope. 
Also, the AGN subsample has a similar slope to 
the quasar sample, which is what one would ex- 
pect if the AGN are obscured quasars. 

4.3. Velocity Dispersion 

The central velocity dispersions of galaxies 
have been found to correlate tightly with the in- 
ferred masses of their central black holes (e.; 



lerred masses ot their central black holes (e.g . , 
Ferrarese fc Merrittl liooot iGebhardt et a l. 2000). 



The latter quantity also scales with the mass of 
the stellar bulge. There is some evidence that the 
most massive bulge systems are more likely to be 
radio loud compared to less massive systems (e.g., 
Laor 2000l : Lacv et al. 2001 ). but in our cases we 
derive median radio luminosities well below what 
one would consider radio-loud objects. As such, 
our analysis pertains more to ordinary systems 
covering a range of AGN indicators (see §[3]). 

The results are shown in Fig. [7j It is clear 
that there is a trend for the most massive sys- 
tems (highest values of a) to have the highest ra- 
dio luminosities. This should not come as a sur- 
prise, based on the correlations between black hole 
mass and bulge mass. As seen in § 14.11 the me- 
dian radio luminosity is apparently not very much 
affected by orientation effects (since the radio lu- 
minosities are comparable for quasars and AGN 
to within a factor of two at the same redshift), 
and provides a rather clean measure of the AGN 
output (at least for the "pure" AGN subsample). 
This then would presumably scale nicely with the 
mass of the central black hole, giving rise to the 
tight correlation between median radio luminosity 




60 70 80 90100 200 
Velocity Dispersion a (km/s) 

Fig. 7. — Median radio luminosity as function of 
central velocity dispersion a of the host galaxy, 
for the AGN dominated (I, pentagons) and star- 
formation dominated (II, stars) subsamples. The 
spread in median redshift is indicated by the color 
coding. The close correlation between the me- 
dian values of the velocity dispersion and the ra- 
dio luminosity is apparent, especially for the star- 
formation sample. The small symbols indicate 
where the colored dots would lie if we use its me- 
dian redshift and the redshift - radio luminosity 
correlation of Fig. [5] Note that the offsets are 
along the y-axis. Low a galaxies are underlumi- 
nous in the radio compared to the average redshift 
trend, and the most massive galaxies are over lu- 
minous in the radio (this holds true for both sub- 
samples) . 

and velocity dispersion (see Fig. [7j pentagon sym- 
bols). The least-squares fit to the star- formation 
subsample (star-shaped symbols) is given by: 



£r\SF 



= (2.50 ±0.05) log (er) + (23.61 ±0.10). (2) 

The AGN-dominated sample does not allow for a 
straight line (in log- log) fit, which may indicate 
other dependencies beyond a simple velocity dis- 
persion / galaxy mass scaling. 

One thing we need to check though, is the effect 
the small range in redshift may have on radio lumi- 
nosity. Based on their tight correlation in Fig. [5j 
a small change in median redshift may account for 
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most of the trend we see in Fig. [7l To this end, we 
calculate the median redshift of the sources used 
for each velocity dispersion bin, and applied the 
redshift - radio luminosity fits from Fig. [5] to infer 
the radio luminosity that way. These values are 
overplotted as the smaller black symbols (so each 
a bin has two radio luminosities associated with 
it). Now, if the large and small symbols would 
closely track each other, then the correlation be- 
tween velocity dispersion and radio luminosity is 
spurious due the sample being flux-limited. How- 
ever, since this is not the case, the radio luminosity 
- velocity dispersion correlation is real. 

Note that the offset in radio luminosity between 
the star-formation dominated and "pure" AGN 
subsamples is even larger than in Fig. [5l at about 
an order of magnitude compared to a factor of 
about 3 (for z < 0.1). Clearly for some objects at 
similar redshifts and similar mass (velocity disper- 
sion), star- formation can account for 90% of the 
total radio luminosity. This would also suggest 
that the correlation for the star-formation subsam- 
ple is driven by the bulge mass - a relation (higher 
a equals a more massive galaxy, which translates 
into a higher absolute amount of star-formation 
for a given D n (4000) value). The radio luminos- 
ity for the "pure" AGN subsample (pentagons in 
Fig. [7|) on the other hand, scales with the mass 
of the central black hole, which itself scales with 
the bulge mass (which is directly related to the 
velocity dispersion). 



Radio Luminosity log(ergs/s) 
36 37 38 39 40 41 42 
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27 28 29 30 31 32 33 
Radio Luminosity log(ergs/s/Hz) 

Fig. 8. — Comparison of the correlation between 
[O III] and radio luminosities (see Fig.0 note that 
the axes are switched) to values from the litera- 
ture. The sample of Best et al. (2005) contains 
SDSS objects with a radio flux of > 5mJy and 
are represented by the small open squares (color 
coded with redshift). The offsets in radio lumi- 
nosity between our median values and theirs can 
exceed 3 orders of magnitude. It should be noted 
that our sample of 80,156 includes the Best et al. 
sample, underlining the severe bias one introduces 
by selecting the b rightest radio sour ces. The solid 
line is the fit from I Wills et al.l |2004) on local low- 
luminosity AGN. 



5. Discussion and Summary 

We note the correlations of the median values 
of the [O III] and radio luminosities in § 14.21 and 
Fig. [SI in particular the similar ratios for the AGN, 
star-formation, and quasar subsamples. Given 
that the mechanisms responsible for the radio 
emission are thought to be different between the 
AGN and star-formation subsamples, this comes 
as a surprise. Our work is not the fir st that that 
touch es on this subject. For instance, IWills et al.1 
(|2004l ) worked on a sample of nearby (z < 0.05), 
low-luminosity AGN taken from a 2Jy radio sam- 
ple. They found an [O III] - radio luminosity slope 
of 0.85, based on a sample of a few dozen source^- 
This value is comparable to the AGN value (0.87, 



3 We converted their 5GHz radio flux densities to 1.4GHz 
values using a constant spectral radio index of —0.5 



see Eqn. 1), but that may be coincidental since 
their sample is so much brighter than ours (com- 
pare the solid line in Fig. [5] to our median values); 
even though the Wills et al. sample is consid- 
ered low radio luminosity, it is still more than 3 
orders of magnitude brighter than our median val- 
ues. They are probing the brightest part of the ra- 
dio source population. This can be put into more 
p erspective by com paring our sample to the one 
of lBest et al. (2005). Their paper is also based on 
FIRST and SDSS data. Actually their data-set is 
wholly contained within ours. 

In Fig. [SJ we replot the data from Fig.[S]onto the 
data from Best et al. We only include sources from 
their sample (open square symbols) that are both 
in our sample and have detected [O III] emission 
(Fig. 9, lower right panel from Best et al. includes 
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a set of [O III] non-detections). Our median [O III] 
values cover the exact same range as the Best et al. 
sample; the main difference is the radio luminosity. 
Since our data represent median values, it implies 
that a full 50% of the sources are located to the left 
of our data-points. The fact that the Best et al. 
data can be up to more than 3 orders of magnitude 
brighter (in the radio) illustrates the severe radio 
selection bias of their sample. It also makes it clear 
that one has to be careful in interpreting possible 
correlations between [O III] and radio luminosities. 
The Best et al. sample would result in a much 
steeper relation. 

Since our sample has not been selected on ra- 
dio properties, and since we include all sources 
in the median calculations, we consider our me- 
dian [O III] - radio luminosity correlation more 
robust than previous work based on radio-selected 
sampled- Consequently, one can use the median 
radio luminosity for the star-forma tion subsample 



to estimate the star-formation rate. Hopkins et al 



([20031 ) compares various star-formation indicators 
using SDSS data, and incorporates the radio- 
derived star-formation rate method of lBelll (|2003r ). 
The inferred star-formation rates are listed in Ta- 
ble 3, and range from less than 1 to about 20 
Mq per year. These modest star-formation rates 
are consistent with values derived through, for in- 
stance, Ha luminosities of star-forming galaxies 



that have been detected by FIRST (Hop kins et al 



2003). The listed values should not be taken at 
face value, however. On the one hand there is 
an unknown AGN contribution to the overall ra- 
dio luminosity which inflates the star-formation 
rate, and on the other hand our median stack- 
ing technique is insensitive to extended (beyond 
~ 6 arcseconds) radio emission, causing it to un- 
derestimate the actual radio luminosity and star- 
formation rate. 

In fact, we measure slightly different sizes for 
the AGN and star-formation subsamples: 5.9" ± 
0.3" by 6.4" ± 0.4" and 5.2" ± 0.2" by 5.7" ± 0.2", 
respectively. The AGN sizes are comparable to the 
median size for SDSS quasars of 6.4" by 7.0" (also 
along RA and DEC), as given in Paper I. This 
implies that there is no extended radio emission 



7 That is not to say that our sample is without biases: galax- 
ies in our sample have to have measurable emission-line 
ratios. 



detected for the star-formation subsample, while 
some has been detected for the AGN subsample 
beyond the ~ 5.4" FWHM of the FIRST point- 
spread- function. 

5.1. Summary 

Combining a large number of objects across a 
range of AGN activity from the SDSS with the 
FIRST radio survey allows us to correlate vari- 
ous sample statistics with radio emission at lev- 
els that are too low for individual sources to be 
detected. We are able to measure median radio 
emission down to a few tens of /xJy, well below the 
FIRST detection threshold. 

We find the following quantities to correlate 
strongly with the median radio luminosity: 

1. D n (4000) index / star- formation rate. This 
index traces the amount of star-formation in 
a given galaxy by measuring the amount of 
blue light relative to the red. Low values 
of the D n (4000) index correspond to more 
active star-formation, and are found to have 
the highest median radio luminosities (more 
than an order of magnitude larger than high- 
index objects). 

2. [O III] luminosity. Higher levels of [O III] 
luminosity correspond to higher median ra- 
dio luminosities. This is true over at least 
3 orders of magnitude. The [O III] - radio 
luminosity correlation for the "pure- AGN" 
subsample is consistent with the same rela- 
tion for quasars, but extended towards lower 
luminosities. The highest levels of [O III] 
and median radio emission are seen among 
the "pure-star- formation" subsample, with 
levels comparable to the quasar ones. 

3. Central velocity dispersion a. Generally, 
higher values of a (i.e., larger galaxies) cor- 
respond to larger median radio luminosi- 
ties. However, there are clear differences 
between the AGN and star-formation dom- 
inated subsamples. For a given median ra- 
dio luminosity, star-forming systems have a 
much smaller a than AGN-dominated sys- 
tems. Conversely, for a given galaxy size and 
a, star-formation dominated systems can be 
an order of magnitude more luminous in the 
radio than the AGN-dominated systems. 
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We furthermore find that the AGN contribution 
to the overall median radio luminosity is roughly 
constant in these galaxies at the 5 x 10 27 ergs s _1 
Hz -1 level. The other component is very tightly 
correlated to the star-formation rate (as derived 
from the optical spectrum). This also implies that 
if one is considering a system without an AGN, 
its star-formation rate can be accurately assessed 
using the radio luminosity of the system. There is 
no indication that this correlation cannot be ex- 
tended to star-formation rates below 0.1 M© / yr. 

It is clear that by studying (median) statistical 
properties of carefully selected and representative 
samples, one can infer subtle trends that other- 
wise are either too small to detect, or are masked 
by source intrinsic variations. Our median radio 
luminosity measurements are at least a few orders 
of magnitude below typical studies (see Fig. [5J, 
and provide a much less-biased view of the radio 
properties of galaxies. Methods such as that used 
in this paper hold great potential when applied to 
future, large-scale multi-wavelength surveys. 
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Table 1 

Median Results for Star-formation Binning 



^median 


SFR 


1<J 6 1, -'-'radio ^ 


F i 

± peak 


D 1 40001 




[AT / , rT ,l 

L M © / yrj 


[ergs/s/Hz] 






U.U / 


n 1 5 


z / .oyinu.uz 


01. ( zno.u 




0.0845 


0.47 


27.90±0.02 


80.4±3.6 


1.754 


0.0819 


0.79 


27.96±0.02 


91.0±3.8 


1.728 


0.0832 


1.12 


28.06±0.01 


101.8±3.9 


1.704 


0.0842 


1.47 


28.14±0.01 


114.2±3.4 


1.697 


0.0865 


1.85 


28.24±0.01 


128.0±3.2 


1.683 


0.0908 


2.26 


28.26±0.01 


127.1±3.8 


1.671 


0.0946 


2.72 


28.36±0.01 


139.3±3.6 


1.653 


0.0994 


3.26 


28.41±0.01 


143.9±3.5 


1.639 


0.1045 


3.90 


28.54±0.01 


172.9±3.6 


1.625 


0.1081 


4.68 


28.62±0.01 


190.7±3.8 


1.612 


0.1151 


5.65 


28.71±0.01 


193.5±3.6 


1.595 


0.1193 


6.97 


28.80±0.01 


224.6±3.4 


1.584 


0.1292 


9.08 


28.92±0.01 


242.9±3.5 


1.565 


0.1481 


14.1 


29.15±0.01 


316.8±3.9 


1.527 



Note. — Both the radio luminosities and flux densities have 
been corrected for the snapshot-bias (see text). The star- 
formation rate is taken from Brinchmann et al. (2004) 



Table 2 
Results for D n (4000) Binning 



D n (4000) 


•^median 


^mean 


log(_L ra dio,flux) 


log(L ra dio,lum) 


log(£[om]/£©) 


log(£[oii]/£©) 








[ergs/s/Hz] 


[ergs/s/Hz] 


2.044 


0.0724 


0.0818 


28.15 


28.05 


6.08 


6.50 


1.942 


0.0826 


0.0927 


28.21 


28.10 


6.12 


6.54 


1.885 


0.0873 


0.0963 


28.20 


28.05 


6.18 


6.57 


1.836 


0.0928 


0.0994 


28.26 


28.10 


6.28 


6.64 


1.791 


0.0951 


0.1009 


28.24 


28.06 


6.30 


6.71 


1.747 


0.0960 


0.1027 


28.30 


28.14 


6.39 


6.79 


1.705 


0.0972 


0.1034 


28.36 


28.21 


6.53 


6.91 


1.665 


0.1008 


0.1063 


28.42 


28.28 


6.65 


7.01 


1.624 


0.1014 


0.1077 


28.48 


28.35 


6.79 


7.17 


1.584 


0.1031 


0.1079 


28.56 


28.44 


6.92 


7.30 


1.541 


0.1053 


0.1095 


28.66 


28.55 


7.04 


7.39 


1.496 


0.1069 


0.1106 


28.73 


28.62 


7.17 


7.53 


1.446 


0.1114 


0.1147 


28.86 


28.79 


7.34 


7.69 


1.384 


0.1143 


0.1178 


28.99 


28.94 


7.54 


7.86 


1.276 


0.1201 


0.1244 


29.21 


29.19 


7.85 


8.11 



Note. — log(L ra dio : flux) and log(L ra dio,ium) are luminosities based on flux-stacking and luminosity- 
stacking, respectively. The bolometric luminosity of the sun (L©) is 3.82 6 x 10 33 ergs s _1 . 
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Table 3 

Inferred Star-formation Rates - Star-formation subsample 



-^median 


log (X ra dio) 


SFR 




[ergs/s/Hz] 


[M / yr] 


0.075 


27.60 


0.45 


0.095 


28.10 


1.07 


0.100 


28.34 


1.61 


0.101 


28.46 


1.97 


0.106 


28.57 


2.37 


0.108 


28.72 


3.06 


0.112 


28.78 


3.38 


0.115 


28.91 


4.49 


0.119 


28.99 


5.40 


0.123 


29.08 


6.64 


0.128 


29.23 


9.38 


0.134 


29.34 


12.1 


0.149 


29.65 


24.7 



Note. — For reference, an L* galaxy 
corresponds to log (£ ra dio) ~ 28.81 and 
an SFR of 3.57 M / year. 
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